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We study the occurrence of 50(5) symmetry in tlie low-energy sector of two-cliain Hubbard-lilce 
systems by analyzing the flow of the running couplings (g-ology) under renormalization group in 
the weak-interaction limit. It is shown that SO{5) is asymptotically restored for low energies for 
rather general parameters of the bare Hamiltonian. This holds also with inclusion of a next-nearest- 
neighbor hopping which explicitly breaks particle-hole symmetry provided one accounts for a dif- 
ferent single-particle weight for the quasiparticles of the two bands of the system. The physical 
significance of this renormalized 5*0(5) symmetry is discussed. 
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Recently, it was shown that both the two-dimensional 
(2-D) Hubbard and the t~J model enjoy an approximate 
50(5) symmetry |l|,y], unifying antifcrromagnctism (AF) 
with d-wave superconductivity (dSC) g]. This symme- 
try principle gives a definite microscopic description of 
the AF — > dSC transition as the chemical potential is 
varied. From the 5*0(5) multiplet structure verified by 
exact cluster diagonalization, one can see how the 50(5) 
superspin vector is rotated from AF to dSC direction and 
show that at the critical chemical potential, the energy 
barrier AE between AF and dSC states is an order of 
magnitude smaller (~ J/10) than the exchange coupling 
J, i. e., the natural parameter in the model. This finding 
is clearly of importance: while it is well-established that 
both t — J and Hubbard models reproduce very success- 
fully the "high" and "medium" -energy physics of order 
U and J ~ t^ /U of the cuprates, the low-energy content 
of order of the superconducting gap has so far eluded 
theoretical investigations. The variance AE* of the mul- 
tiplet splitting is a well-defined measure of "how-good" 
the 50(5) symmetry is realized in the bare model con- 
sidered for a given system size. It seems then a natural 
question to ask, whether this deviation from exact 50(5) 
symmetry remains small or even vanishes if one goes to 
the infinite-volume limit and to lower energies, i. e. un- 
der renormalization-group (RG) flow. 

In this Letter, we demonstrate that the low-energy 
regime of rather general Hubbard-type models including 
finite-ranged interactions (provided they are weak) be- 
tween one and two dimensions, i. e. ladders ^ , is indeed 
dominated by the scaling towards an 50(5) -invariant 
model. This is a remarkable result, since it is the first 
model, which is non-50(5) invariant at the bare starting 
(microscopic) level, where the existence of 50(5) sym- 
metry is proven for low energies g]- We first show that 
this holds for the particle-hole- (ph) symmetric case with 
nearest- neighbor hopping only n. In addition, we con- 
sider the effect of a next-nearest-neighbor (intrachain) 
hopping i2 which produces an explicit breaking of the 
ph symmetry. We demonstrate that also in this case the 
system becomes 50(5) symmetric for low energies (at 
least up to order t\) provided one considers a renormal- 
ized S0{5) transformation M, which takes into account 



a different renormalization of the single-particle weight 
of the bonding and antibonding bands. This result is of 
importance because of two points: (i) it sheds light on 
the effect of the longer-range hoppings in general on the 
fate of 50(5) symmetry. This issue is also under inten- 
sive discussion in the case of 2-D systems |@,||9| : ^2 is 
known to strongly affect AF correlations and the Fermi- 
surface topology in the cuprates. (ii) In addition, the 
renormalized 50(5) representation introduced here for 
the first time, is likely to be realized in a significantly 
larger class of physical systems allowing, for example, 
for asymmetries in the AF- and dSC- phases, such as 
different transition temperatures. 

Specifically, we consider two coupled chains in the band 
representation with total "low-energy" action 5 = 5o + 
5/, where the non-interacting part Sq can be written at 
a given point r in the RG fiow as 
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Here, Oko- (C'^o-) ^'"^ Grassmann variables associated with 
the destruction (creation) of a fermion, and a is the spin. 
k = {iw, /cx, fell , z^} is a shorthand notation for the Mat- 
subara frequency ioj, and the momentum perpendicular 
(fcx = (0,7r)) and parallel (fc||) to the chain direction. 
The latter momentum is measured relative to the Fermi 
point V kp_k± associated with the "band" k± with Fermi 
velocity Vk^n and v = ±1 refers to right- and left-moving 
fermions, respectively. This action is restricted to modes 
with |fc|| I < A with A — Aqb^^. This weak-coupling RG 
method has previously been applied to obtain the phase 
diagram of the two-chain Hubbard model [J] . In order 
to study the occurrence of 50(5) symmetry in Hubbard- 
like models, it is convenient to rewrite the interaction 
part 5/ of the action in terms of 50(5) -invariant and 
50(5) -breaking terms Q . Defining the 50(5) spinor 
as in Ref. [fll 
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where k stands for {— iw, —k±-\-TT, —k\\, v}), the interact- 
ing action 5/ can be shown to be expressable in terms of 
the 4x4 charge- rotation Dirac matrix F^'"^ pdl]: 



Sr- 



2(3N 



E '5o(---)<(l + «(---r')*k. 



X ^l^ (1 + b{- ■ oris) *k4 + [kuh) ^ {k2,h) ■ (3) 

Here, (• • •) represent the sets of variables on which 
the couplings go, a, and b depend. As usual, each 
coupling can be considered as dependent only on the 
Fermi momenta closest to where the corresponding pro- 
cess takes place, i. e. the (• • •) are labeled by 
{h'i,k±i;v2,k±2W3,k±^;h'4,k±4). Moreover,^ denotes 
a sum with conservation of frequency and lattice momen- 
tum. 

The S0{5) -symmetric part Sj of Sj is given by 
Eq. d) with a{- ••) = &(•••) = |l|,|ll. For a gen- 
eral S0{5) -invariant action S} it can be shown that 
one can restrict oneself to the seven independent cou- 
plings g^^\ g'i'°~'"\ g'i\ 90 '^"^K go\ 3o"'°~''^ and 
5q ' ^ , defined in analogy with the g-ology formalism 
0-|l6|. The 5*0(5) -breaking term SJ^^ of the interacting 
action with ph symmetry can be demonstrated to be the 
term proportional to go{- ■ ■)o(- ■ ■)b{- ■ •) in (j^) and thus 
we define the corresponding 50(5) -breaking couplings 
as gi{- ■ •) = go{- ■ •) a(- • •) &(• • •). In this case, we can re- 
strict ourselves to only 5 independent couplings (it can be 
shown that the others are redundant), namely, gl , gi , 
^(1)^ ^(it:0-7r)^ ^^^ gO--^o~^) g^ At half-filling and with 
^2 = 0, the Hamiltonian is ph symmetric, since the veloc- 
ities of the two bands Vfcj^^o a-nd Vk^=Tr are equal. Since 
the ph-breaking terms in (0) are proportional to Fis, one 
can set in the ph-symmetric case a(- • •) = —b{- ■ •), and 
consider the RG flow of the couplings ^q and gl only. 

To begin with, we have evaluated the RG equations 
for the (7q and g^ couplings |p6| at one loop by using 



the standard 5-ology procedure (cf. Refs. [|l5|J4|j, m- 
cluding the interband umklapp processes. As already 
shown for the two-chain case |4| , the system always flows 
to strong coupling, i.e. the g's diverge at a value of 
T = Tc (X l/.g(r = 0), g{T) being the scale of the interac- 
tion (proportional to the maximum of all g^ (t)). This 
signals an instability towards some gapped state. Nev- 
ertheless, the striking new result is that even in a non- 
50(5) -invariant system, like, e.g., the Hubbard model, 
the 50(5) -invariant couplings ^q dominate with respect 
to the symmetry-breaking couplings gl , when approach- 
ing Tc [17| . This can be seen from the ratio of the max- 
ima of these two types of couplings, gi°"^{T)/g^°'^'{T) 
going to zero, as shown in the inset of Fig. |l|. Here, 
5j (r) is defined as the largest absolute value, and 
thus the scale of the couplings of a given type i {i = 
0,1, ph) at a given r. This result implies that the low- 
energy modes of the system can be described by an ef- 
fective SO{5) -sym,m,etric action, at least for sufficiently 
small g{T = 0) fl^fl. In fact, we have verified that this 
occurs for very general values of the Hamiltonian, includ- 
ing longer-ranged interactions. 

A next-nearest-neighbor hopping ^2 breaks ph symme- 



try explicitly and requires the introduction of a ph break- 
ing interaction Sj^'. Here, a(- • •) 7^ — &(• • •) in Eq. (||), 
and thus we need extra couplings gphi' ■ ■), which we have 
defined as gph{- ■ •) = 5o(- ' O («(■ • •) + K ' O) /2- In this 
case, one can show that the couplings can be restricted 
to ffit\ gf^, and gi]^ mm. The initial (r = 0) source 



Jph ' yph ' "^ ^ Uph 

of ph-symmetry breaking for ^2 7^ stems from the non- 
interacting part of the action Sq, due to the difference of 
the Fermi velocities AVb of the two bands. In the follow- 
ing, we will show that 50(5) symmetry is restored (at 
least up to 0(^2)^), at low energies, even in the presence 
of this ph- (and thus 50(5) -)breaking term. 
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FIG. 1. 

RG flow of the 50(5) -breaking terms 3p$™'"'V5o""'"'' (fuU), 
3<7pr''"'/5o""'"' (dashed), AK/AK=o (dotted), and AZ^' 
(dash-dotted), as function of r = — log(A/Ao) for [7 = 1, 
t = t^ = \,t2 = 0.5 and half fiUing. Here, AK- = Vot — Virr, 
and AZ-i = {{ZZ,^IZ-^f - l) / (T/o.=o/Kr=o - 1). The 
inset shows Sgj ""^ /(?g "^ vs tU for ±2 — 0. 

These results are obtained on the basis of two comple- 
mentary RG calculations. Calculation (i) considers the 
RG flow of the self-energy parameters Vk^r and Z'^^^ 

at two loops, and of the coupling parameters ffj (t) at 
one loop, taking the t dependence of all the parameters 
at each RG step fully into account. This first calcula- 
tion (i), although not rigorously controlled (see below), 
is motivated by the fact that we are interested in study- 
ing the RG flow of the self-energy, which is the leading 
symmetry-breaking term when ^2 is included pq] . In a 
second calculation (ii), we will show how our main results 
about the renormalized 50(5) symmetry obtained within 
this first procedure can be achieved also in an alternative, 
more controlled way, where we consider only the renor- 
malization of the g^ . Nevertheless, the first calculation 
(i) is instructive, in order to provide a physical interpre- 
tation for the single-particle renormalization factors Z as 
discussed in the conclusions. Indeed, in procedure (i) the 
Z factors, and thus the renormalized 5*0(5) transforma- 
tion, derive naturally from the RG flow, while in (ii) they 
are introduced right at the outset. 

In calculation (i) , the relevant part of the renormalized 
action has the form of Eq. (nl) with r-dependent Fermi 
velocities and single-particle weights. The flow of these 
parameters is shown in Fig. ^ As the bare (t = 0) 



Hamiltonian, we take the half-filled Hubbard ladder with 
isotropic intrachain and interchain hoppings t = t± = 1, 
next-nearest-neighbor hopping t2 = 0.5 (corresponding 
to Ay^=o « 1.9), and [/ = 1 ||l|]. Actually, we have 
verified that the results we are discussing below are rather 
general and hold also in the presence of anisotropy tj_ ^ t 
and nearest-neighbor interactions (^ —U). The t2 = 
case P], discussed above, is plotted in the inset for 
comparison. 

For the ^2 7^ case, AVr (dotted hne) flows to zero, but 
AZ~^ (dash-dotted line, initially zero) scales to unity. 
Therefore, the initial asymmetry between the bands due 
to the different Fermi velocities is transferred into a dif- 
ference in the single-particle weights Z, such that for 
large r Zq^/Z:^^ — > \/VoT=o/K-r=o- In order to restore 
the coefficient of the ioj term in (|l|) to unity, the stan- 
dard procedure |15|| is to reabsorb this renormalization 
into the definition of new Grassmann variables Cko- and 



to set y Zj. ^Ckcr — C^a- This standard procedure is dic- 
tated by the requirement to identify the canonical Fermi 
operators with correct anticommutation relations, as will 
be discussed at the end. 









2 


4 T 


6 


8 


10 


1.0 














RG 














0.5 




(b) 








A 














_ _ _ _ 


\ 



1.0 



tU 



FIG. 2. SO{5) -breaking couplings 33l""'"'V5o'"''''' (fuU) 
and Sg^"^^ / Qq"^^ (dashed) as a function of r. (a) shows 
the results of the RG procedure (i) and (b) shows the results 
of (ii) as defined in the text. The parameters are as in Fig. 1. 

In this way, the non-interacting part of the (renormal- 
ized) action will again be symmetric under exchange of 
the two bands (and thus 15*0(5) symmetric in the new 
fields). This transformation, however, also affects the in- 
teraction part, and one should consistently redefine the 
renormalized 5*0(5) spinor in (0) to ^k, whereby the 
Okcr are again replaced with the Oko-- The couplings de- 
fined in this way are of course different from the origi- 
nal ones and we will distinguish them with a tilde, i.e. 



9. 



-^ gl . The remarkable result is that the transforma- 



tion which makes the non-interacting part of the action 
S0{5) -symmetric also restores 5*0(5) in the interacting 
part. This is demonstrated in Fig. Ha, which plots the 



ratio of the g. 



■(max) 



as a function of the fiow parameter r. 



We note that the non-50 (5) couplings g)^°'^' and g^ 



all flow to zero (relative to the ^g "^ ). Thus, at large r, 
50(5) symmetry is restored for low energies in the "~ " 
basis. However, at the energy scale where AV^ starts 
to decrease and /S.Z~^ starts to become finite (t ~ 7 in 
Fig. P, the renormalized couplings can be shown to be- 
come large and the weak-coupling expansion is no longer 
controlled, as anticipated. 

To support this physically appealing yet uncontrolled 
calculation, we verify, in terms of a controlled RG calcu- 
lation |17| (i. e., at one loop), that 50(5) symmetry is 
indeed recovered at least up to 0{i^). This alternative 
derivation clarifies why the single-particle weights renor- 
malize proportionally to (V^r^o)""'", as obtained asymp- 
totically in the two-loop calculation. In this RG proce- 
dure (ii), we start from the action Sq + Sj and carry 
out the transformation T on the Grassmann variables 
right at the outset, where T is defined as TCiuj^kn,k± — 

l/V"fc±C'i^',fc|i',fc_L ■^it^^ *^' = *^/wfc_L and k\\' = k\\ Uk^, 
and Uk^ = v^Wr [@' ^^^^ a transformation, which is 
always possible with Grassmann variables, is motivated 
by our first calculation (i). By changing the sum over 
iu and fcy into a sum over iuj' and fcy separately for 
each band, the non-interacting part of the action again 
recovers its explicit 50(5) symmetry. Furthermore, by 
defining a new 50(5) spinor ^k in terms of the C, we 
obtain new couplings 5q , gi , and gA, as in step (i). In 
Fig. 0b, we show the corresponding RG fiow of the ra- 
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With increasing RG parameter r the 



ph-symmetry breaking term 5^™"^ /ffg °^ vanishes (full 

line), while the 50(5) -breaking term ^^™"^Y^g '"^' goes 
to a finite but rather small value |2y] (dashed line). The 
50(5) symmetry thus is recovered up to a very high de- 
gree of precision for low energies p]. In contrast with 
the results of procedure (i) , the result of (ii) is controlled 
for small g{T = 0) |l^,^. In this way, we have shown 
in a controlled way that 50(5) is restored for low en- 
ergies at least up to order t| for small g{T — 0). Our 
two-loop calculation (i) further suggests that even this 
50(5) -breaking term of order t^ might be removed by 
the self-energy renorrnalization. 

The renormalized 50(5) symmetry introduced here, 
and the related renormalization of the single-particle 
weights Zk^T, can be understood in terms of a simplified 
scheme, which renormalizes the Hamiltonian, by restrict- 
ing the Hilbert space to a subspace with energy uj smaller 
than a certain cutoff luq oc Ag exp(— r). In the restricted 
subspace, the total integrated spectral density (which we 
identify with Z^^t) will be less than one. Since the spec- 
tral sum rule identifies Zk^r with the anticommutator of 
the Fermi operators Ok, the canonical Fermi operators 
with anticommutator equal to 1 in this subspace are the 
transformed field operators Ok introduced above p3|. 

In conclusion, we have shown that the effective low- 
energy action (or Hamiltonian) of a ladder with weak 
interaction is asymptotically 50(5) symmetric p|. With 
the inclusion of a next-nearest-neighbor hopping ^2 the 



action is invariant under a generalized SO{5) transfor- 
mation 10], which performs a "stretched" 5*0(5) rotation 
of the order parameters. Physically, this SO{5) symme- 
try may be present in the low-energy sector of a larger 
and more generic class of physical systems than the ordi- 
nary 50(5) . Moreover, since this stretched rotation does 
not conserve the norm of the superspin (order-parameter) 
vector [B|, a renormalized 50(5) theory can possibly ad- 
mit asymmetries between the antiferromagnetic and su- 
perconducting phases, like for example the difference in 
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In order not to confuse the reader, we have chosen for 
this new transformation the same symbol " " as for 
the one performed in the procedure (i), since they have 
the same physical meaning. The two transformations dif- 
fer for large r by (a) an overall constant, and (b) by a 
scaling of the frequencies and of the momenta. Notice, 
however, that due to their frequency- and momentum- 
independence, the gl couplings defined in (ii) are equiv- 
alent to the gl defined in (i) apart for a common factor. 
We have verified that, as a function of f2 and for different 
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The only remaining explicitly non-SO(5) invariant term 
in the renormalized action is due to the different mo- 
mentum and frequency conservation for the k± — and 
k± = ir bands arising from the T transformation. How- 
ever, one can show that this term also breaks the SO(5) 
symmetry by an amount of order f^ at most. 
In the case of the ladders, where Zk^r vanishes at Tc due 
to a spin gap, this argument can be applied up to, say 
T — T* for which Zkj^T is still appreciable and still the 
SO(5) symmetry is achieved (see \V^) . 



